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Abstract 

I discuss the slow convergence of Discretized Light Cone Quantization (DLCQ) in the small mass 
limit and suggest a solution. 

1 Introduction 

"Light-front field theory" is an attempt to apply Hamiltonian techniques to relativistic field theory; the 
idea is to quantize a theory on a surface of constant x+ = {x^ + x^)/\/2 and calculate eigenstates of the 
invariant mass squared operator P^^P^. In order to partially regulate a theory, one can simply impose 
(anti-)periodic boundary conditions in the coordinate; this is known as "Discretized Light Cone 
Quantization" (DLCQ) [Q since the momenta conjugate to x~ become discrete. One can use DLCQ to 
construct a basis of states for use in numerical calculations. 

Why would one want to use DLCQ in a large scale numerical calculation? Ideally, one would employ 
a basis of polynomial- type wavefunctions to solve the bound state equation Q|. With such a basis, 
convergence of the eigenvalues is typically exponentially fast. However, the practitioner faces several 
difhculties in applying this approach to a many particle calculation. In general, the resulting basis is 
not orthogonal and the resulting Hamiltonian matrix is no longer very sparse. Second, matrix elements 
become difficult to compute since they involve complicated integrals. It is not clear that the faster 
convergence of such a basis can compensate for these added difficulties. 

In contrast, a DLCQ basis is orthogonal and the corresponding Hamiltonian matrix is sparse. Also, 
matrix elements are easily calculated; no integrals need to be evaluated. This is important since the 
time required to calculate matrix elements is currently the limiting factor in large calculations. Finally, 
as a regulator, DLCQ has an appealing physical interpretation in terms of (anti-)periodic boundary 
conditions. 
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Figure 1: Lowest eigenvalues of the 't Hooft equation as a function of mass /i calculated using a 
polynomial wavefunction basis . 

However, DLCQ has one distinct disadvantage: it converges quite slowly in the small mass limit. 
This is especially of concern for the transverse lattice Q where the continuum limit is precisely the 
small mass limit. In particular, the lowest eigenvalue and measurements of the string tension are strongly 
affected by this slow convergence. As we shall see, the source of this slow convergence comes from the 
end-point behavior of the wavefunction which is handled poorly by DLCQ. Using the 't Hooft equation 
as an example, we will study this behavior and suggest a solution that can be easily applied to other 
theories and large numerical calculations. 

2 The 't Hooft equation 

The 't Hooft equation [|| is the bound state equation for two quarks in the large iV limit of (1+1)- 
dimensional QCD: 

\x i-xj Jo [x-y) 
where M is the Lorentz invariant mass eigenvalue associated with eigenfunction ^i(x). A Cauchy principle 
value prescription is assumed throughout. The spectrum is shown in Fig. ^ Following 't Hooft, let us 
analyze the endpoint behavior of the wavefunction. Assume 

lim ^{x) oc x'^ . (2) 

For small x, the bound state equation becomes 

= V' - - f dy (1 + 0{y)) + O(x^) . (3) 

Jo [x-yr 
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Substituting z = y / we obtain, in the a; — s- limit, 

,2 /■!/ 



^ - 1 + 7r/3cot(7r/3) . (5) 
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Convergence of the integral ^ implies < /3 < 1. For small mass 



/^=^ + • (6) 



In addition, one can derive an expression for the lowest eigenvalue in this limit 



M^^^-^ + 0{^?) . (7) 

3 Analysis and discussion 

Exact solutions of the 't Hooft equation are not known and one must resort to numerical techniques. 
For the purpose of comparison, we will introduce a basis of wavefunctions of the form x^{l — x)^Pn{x) 
where Pn{x) is some polynomial of order n and /3 is given by Eqn. (^. Such a basis provides very 
quick numerical convergence [Q. Later, when I quote "exact eigenvalues," I will use results from such a 
calculation. 

Alternatively, we can use DLCQ to solve the bound state equation. We discretize the momentum 
fraction x as i/K where i G {1/2, 3/2, . . . ,K ~ 1/2} for integer-valued "harmonic resolution" K. Physi- 
cally, this discretization corresponds to applying anti-periodic boundary conditions in the x~ direction 
where ii' ^ cx3 is the continuum limit. The DLCQ version of the bound state equation ([^) is 

Although a proof is lacking, numerical evidence suggests that the DLCQ eigenvalues converge as l/K"^^. 
For example, in Fig. |^, we plot the lowest DLCQ eigenvalue as a function of K together with a numerical 
fit to 

,.2 1.25943 1.15609 1.01505 

= 0.779141 - + ^ + . . . (9) 

where the "exact eigenvalue" (from a polynomial basis calculation) is A'P — 0.77914. 

Note that the DLCQ bound state equation does not make use of our knowledge of the endpoint 
behavior of the wavefunction. This is, in fact, the source of the slow DLCQ convergence. We can include 
the endpoint behavior by adding and subtracting a term from the 't Hooft equation: 
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M'm{x)= ti'^{x)[- + - \+g'^{x)l{x)+g' dy ^-^ (10) 

\x l-xj Jo {x-yY 
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Figure 2: Convergence of the lowest eigenvalue with K for j3 = 0.1 (/i = 0.181981). The dots are from 
DLCQ and the heavy line is the exact result. 



where 



l{x) 



1 1 
dy — 



( ya-y) y 
{x - yY 



(11) 



The DLCQ version of this equation is 



K 



(12) 



where we evaluate the integral I(a:) exactly]^ Based on the endpoint analysis, the integrand in Eqn. ( [l0| ) 
should vanish when both x and y are near one of the endpoints. Thus, the corresponding discretization 
( |l^ ) should have small errors in this region. Wc will refer to this form of the bound state equation as 
"improved DLCQ." 

Let us see how this affects convergence of the spectrum. In Fig. ^, we see that the DLCQ results 
converge slowly and, in fact, have an incorrect functional form for M"^ versus /i in the small ^ limit. On 
the other hand, the improved DLCQ eigenvalues are nearly identical to the exact result. However, this 
absence of finite K corrections is the result of a fortuitous accident: for small /i, the lowest eigenfunction 
is « x^{l — x)^ . Thus, the last term of Eqn. ( [l^ ) is almost zero for all i and j and all finite K 

corrections vanish. 

Examining an excited state, where the wavefunction has a more complicated form, will give us a 
better picture of what happens in the generic case. Ordinarily, DLCQ eigenvalues converge as 
due to the endpoint regions x and a; « 1 plus l/K errors due to the l/{x — yY singularity Thus, 
in Fig. |[ the DLCQ eigenvalues of the second excited state are fitted to powers of plus powers 

^In practice, I solve it numerically. 
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Figure 3: The lowest eigenvalue as a function of /i. The black lines are from DLCQ (K = 10 and 100) 
and the heavy line is the exact result. The improved DLCQ eigenvalue (gray line, K = 10) is almost 
identical to the the exact result. 
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Figure 5: The third lowest eigenvalue as a function of ^. The black lines are from DLCQ {K = 10 and 
100); the gray lines are improved DLCQ {K = 10 and 100); and the heavy line is the exact result. 



of 1/K, 



,.2 3.86429 3.88433 20.7358 19.0972 

= 16.0016 - + + ^ + , (13) 

while the improved DLCQ eigenvalues are fitted to powers of 1/K, 

, 24.7124 14.7355 
Af 2 ^ 16.0193 —- + —TPT- ■ (14) 

The exact result is = 16.0016. It is important to note that improved DLCQ removes the errors 
associated with the endpoint regions but docs not significantly affect the l/K errors. Finally, in Fig. ^ 
we see how convergence behaves as a function of /i for the second excited state. 

So far, we have discussed the two particle bound state equation. The generalization to a DLCQ 
calculation with more particles is straightforward. For example, consider the following terms from the 
Schwinger model Hamiltonian ^: 

J,t !, 1 jt J °° fe-1/2 



fc=l/2 fe=l/2 n=l 



, ^ ^^-iiZLfetfo^dU. (15) 
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In the two particle sector, this Hamiltonian produces the DLCQ version of the 't Hooft equation (||). 
The improved DLCQ version of (p^ ) is 

± + 9' ± JkAhdU^ - E ^3^6tMU (16) 

fc=l/2 fc,; = l/2 k,l,m,n=l/2 ^ ^ 

k^l 
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where 



i=i/2 

We have replaced the second term of (|l^) by a four-point interaction since the improved DLCQ sub- 
traction is sensitive to the momenta of both particles. The improved form will remove the slow 
convergence from a many particle calculation. 

4 Conclusion 

There have been some past attempts to quantify the errors associated with DLCQ. The claim in Ref. |^ 
that the leading error is 1/K is certainly incorrect. In Ref. 0, periodic boundary conditions were applied 
to a somewhat different theory. In that case, a leading error of l/ln(/'ir) seemed to explain the data. 

There are a large number of DLCQ studies which have suffered from this slow convergence problem. 
This includes DLCQ studies of (l+l)-dimcnsional QCD ^, studies of the Schwinger model (as 
noted in Ref. [0), dimensionally reduced theories [0, |l^, and the transverse lattice However, one 
should note that there are many cases where the slow convergence problem is not important: It is only 
of concern when fi is small but nonzero (the DLCQ subtraction is correct when /i = 0) and 1/K errors 
can dominate in a many particle calculation (where the size of K is somewhat limited). 
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